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Recently, individual-based models originally used for biological purposes revealed inter- 
esting insights into processes of the competition of languages. Within this new field of 
population dynamics a model considering sexual populations with ageing is presented. 
The agents are situated on a lattice and each one speaks one of two languages or both. The 
stability and quantitative structure of an interface between two regions, initially speaking 
different languages, is studied. We find that individuals speaking both languages do not 
prefer any of these regions and have a different age structure than individuals speaking 
only one language. 
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1. Introduction 

Physicists investigate biological and sociological systems, applying their tools, bor- 
rowed mainly from statistical physics. In the last years, the evolution of languages 
and competition among them gained much interest of non-linguistics, especially 
researchers working in evolutionary systems similar to biological ones. There have 
been several approaches to reveal the similarities between the evolution of biological 
systems and languages at the end of the nineteenth century Numerical models 
simulating the competition of many languages have given new insights into the be- 
havior of agents, for instance on a latticeEEEEl as we ll as i n to the size distribution 
of languages^. A review of them can be found in ref. El. 

Our mod el is based on the well understood Penna ageing model I3BI on a i a t_ 
tice I 9 I 1Q | H | which provides us a possibility to model a sexually reproducing stable 
population. We simplify the model by defining languages as an integer number, that 
is they are not composed of different words. Thus we are able to avoid changes in 
the same language. The parents pass their language entirely to their offspring. In 
order to stabilize their distribution, languages can be forgotten during lifetime. 
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The study of an interface between originally different regions under different 
parameter sets reveals certain characteristics of the geographical distribution of the 
languages on the lattice as well as of the age structure of the agents for different 
languages. 

This articles is organized as follows: the next section explains the main features 
of the model and tries to justify the parameters we use. We present the results and 
the conclusions in the two following sections. 

2. The model 

This section is separated into two subsections in order to provide the reader a small 
review of the Penna model on a lattice as well as to present our modifications 
adapting the model to a system where languages compete with each other. 

2.1. The sexual Penna model on a square lattice 

Each individual or agent comprises two bit-strings (diploid) of 32 bits that are read 
in parallel. After birth, at every time step a new position of both bit-strings is 
read. A bit equal to 1 corresponds to a harmful allele. All individuals have five 
predefined dominant positions where one harmful allele suffices to represent an 
inherited disease starting to diminish health from that age on which corresponds 
to the bit position. At the other positions two set bits are needed to switch on the 
effect of a disease. As soon as the agent reaches an age at which the current number 
of deleterious mutations exceeds the threshold value T = 3, the agent dies. In order 
to have a stable population, every time step an individual dies with the additional 
probability: V = N(t)/N max where N(t) is the actual population size and N max the 
so called carrying capacity representing a maximum population size. After reaching 
the minimum reproduction age R = 10, a female agent searches every time step 
for a male, of age equal or greater than the minimum reproduction age, among the 
central site and its four nearest neighbors to generate two offspring. It selects a male 
on the central site with a probability of 25%, and if it fails, it searches among the 
nearest neighbor sites, at each one with a probability of 25%. The two bit-strings 
of the offspring are built by a random crossing and recombination of the parents 
bit-strings (see ref.^j). Each new bit-string suffers a deleterious mutation (from zero 
to one) at a random position. If the selected bit is already 1, it remains 1 in the 
offspring bit-string. The offspring is placed on a nearest neighbor site of the mother 
even if the site is already occupied (which is different from the usual versions of the 
Penna model on a lattice) . Every time step an agent moves to a randomly selected 
nearest neighbor site with probability p m , if this site is less or equally populated. 
The bit-strings are initialized randomly with zeros and ones at the first time step. 

For a more detailed description of the Penna model and its implementation on 
a square lattice we refer to refs. \2i2i 
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2.2. Competition of language 

For simplicity we define a language by an integer number and not by a bi t-str ing 
which would describe, for instance, different words or an alphabet as in ref. EEI 

Every agent speaks a language I, an individual variable which can have three 
values: I = 1 and I = 2 mean that it speaks language 1 or 2, respectively. The 
third possibility, I — 3, describes the case where an agent speaks both languages. 
Our model contains two parameters dealing directly with these language values. At 
birth an offspring learns the language I of its parents if they speak the same one(s). 
In the case of parents with different values of I the offspring speaks both languages 
(I = 3) with probability pb, otherwise it speaks only one language, each with the 
same probability. The other parameter is the probability pf, for which an agent 
may forget an already learned language. Every time step an agent, which speaks 
both languages, counts the number of surrounding agents speaking language I in 
its neighborhood. This neighborhood is defined by a square of a distance of d sites 
from the central site, for instance the 8 nearest neighbors with d = 1 or 24 with 
d = 2. The central site is not counted. If and only if there is a majority of people in 
the neighborhood speaking language 1 or 2 the agent forgets language 2 or 1 with 
probability pf, respectively. Thus it speaks only the language which dominates in 
its surrounding. The lattice has free boundary conditions. 

3. Results 

In our simulations we restrict ourselves to the following initial conditions: half of the 
lattice of size L x L is filled up with agents speaking language I = 1 and the other 
half with language 1 = 2. We study stability and shape of the interface between 
these two regions. Simulations with randomly distributed languages do not present 
a stable interface: The whole population speaks only one and the same language 
after a few time steps. 

The initial population consists of 10, 000 males and 10, 000 females randomly 
distributed over the lattice with the values of I as described above. The carrying 
capacity is N max = 1, 000, 000 on a 20 x 20 square lattice. 

The simulations show that the interface is neither stable for pb values smaller 
than one nor for low occupation (less than 100 agents per site), which is controlled 
by the carrying capacity P max and the lattice size. After a short time the number of 
agents speaking I begins to fluctuate strongly and finally converges into the stable 
state where only one language is spoken. In general the stability of the interface 
depends crucially on the initial state as well as on the random seed. For instance, 
Figure H shows a simulation with pf = 0.1, pb = 1, p m — and d = 1 presenting 
such an instability. 

We concentrate now on the results of the simulations for different values of the 
parameter pf, leading an agent to forget one of its languages. We fix the other 
parameters to pb = 1, d = 1 and p m = 0. As a function of age, Figure El shows the 
population of bilinguals (agents speaking both languages) divided by the monolin- 
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Fig. 1. Unstable interface: Only agents with I = 1 remain. 



guals (I = 1 or 2) for different values of pf. The function decreases exponentially due 
to the effect that every time step a fraction of the bilinguals becomes monolingual. 
With increasing pf a larger fraction of older bilingual agents forgets one of their 
two languages since in their environment they do not need both. The higher is the 
probability to forget a language, the smaller is the number of older agents speaking 
both languages and thus less offspring with I = 3 are created. Figure depicts the 
mean value of monolinguals and bilinguals during one simulation for different values 
of pf. The number of bilinguals decreases with increasing pf as expected. It seems 
that the fraction of bilinguals decays roughly as power law with exponent —1 for 
higher values of p/ . 

Figure^shows the number of agents with certain value of I versus their position 
in direction x perpendicular to the interface. The number is averaged over the 
direction parallel to the interface and is measured after 10, 000 time steps. We 
observe a quite stable interface between the two regions, each one with one of the 
two languages in majority. The number of agents speaking I = 1 and I = 2 decays 
exponentially at the interface as also reported in ref . El. Interestingly, the number of 
bilinguals is constant over the whole lattice. The shape is not altered by changing 

Pf- 

We increased the number of agents by setting N max = 10, 000, 000 and the initial 
populations to 100, 000 females and 100, 000 males for pt = 1: Now the exponential 
decay at the interface is observed clearly, as seen in Figure 
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Fig. 2. The number of agents speaking language I = 3 divided by the number of agents with 
I = 1, as a function of ages. The number of agents speaking both languages decreases drastically 
with age for large values of p j . 



In our simulations we have also changed the parameter d, defining the number 
of neighbors an agent with I = 3 examines, in order to know which language it can 
forget. The results for large d are the same as for d = 1. 

The distribution of speakers on the lattice for different movement rates p rn is 
shown in Figure We set pf = 0.2 and d = 1. Higher movement rates lead to a 
smoother interface. Thus the exponential decay is weaker for large p m . At very high 
movement rates the interface becomes unstable. 



4. Conclusions 

We present simulations where the population of speakers of two different languages 
are of similar size for at least 10, 000 time steps. This meta-stable state is obtained 
only for a large number of agents per site and initialization of the lattice by dis- 
tributing the speakers of different languages on different halfs of the lattice. We can 
interpret that an interface of speakers in high populated areas, for instance at the 
Canal Street of New Orleans where on one side French and on the other side English 
is spoken, is more stable than in low populated areas. Different languages cannot 
survive for long times if their speakers are not geographically separated. Another 
result of the model is that the fraction of bilinguals relative to monolinguals de- 
creases exponentially with age. Older people, living for long time at the same place, 
do not need a second language. The mean value of the number of bilinguals versus 
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Fig. 3. The mean value of the number of bilinguals seems to decay as a power law. 



the parameter pr to forget a second language shows a power law with exponent — 1. 
\X\(i\ 

The results of ref. are well reproduced in our model: At the interface the 
number of monolinguals decreases exponentially. Surprisingly, the number of bilin- 
guals distributes rather homogeneously over the whole grid. The exponential decay 
of the number of monolinguals becomes steeper for smaller movement rates but 
is left unaltered by the lattice size. Higher movement rates lead to a more homo- 
geneous distribution and can break the meta-stability of the interface. Nowadays, 
globalization gives us the possibility to travel frequently over long distances and to 
stay larger periods at different places on Earth, one of the reasons why languages 
are becoming extinct. 

We presented here the first model for language competition including ageing and 
sexual reproduction and reproduced well the results of other models although they 
are quite different. Numerical agent-based models on the computer yield interesting 
results despite their simplicity, and we think that there will be much more to be 
done in future. 

Acknowledgements 

I am funded by the DAAD (Deutscher Akademischer Austauschdienst) and thank 
D. Stauffer and S. Moss de Oliveira for very useful comments on my manuscript. 



9, 2008 2:41 WSPC/INSTRUCTION FILE Cultgrid 



Simulation for competition of languages with an ageing sexual population 7 




Fig. 4. The frequency of monolinguals saturates for high pf values. The dashed line corresponds 
to I = 1 and the solid line to I = 2. For all values of pf we find an exponential behavior at the 
interface. 
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Fig. 5. Frequency for all I on the interface. Bilinguals arc homogeneously distributed. The number 
of monolinguals decays exponentially at the interface. 




Fig. 6. Frequency of monolinguals on the lattice for different movement rates p m with d = 1 and 
Pf = 0.2. The larger p m is, the smoother becomes the steepness of the exponential decay. 



